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tJ" ' Abstract 

o 

en ■ Matrix elements of the form < 0\Tr g^GG\G > are calculated using the 

lattice QCD Monte Carlo method. Here, |G > is a glueball state with quantum 
I numbers 0^^, 2++, ^ and G is the gluon field strength operator. The matrix 

elements are obtained from the hybrid correlation functions of the fuzzy and 
plaquette operators performed on the 12^ and 14'* lattices at /3 = 5.9 and 
r~| ! 5.96 respectively. These matrix elements are compared with those from the 

QCD sum rules and the tensor meson dominance model. They are the non- 
perturbative matrix elements needed in the calculation of the partial widths of 
J/^' radiative decays into glueballs. 



X 



The glueballs predicted from Quantum Chromodynamics (QCD) are the most ex- 
otic particles from the point of view of the naive quark model, but they are crucial to 
the QCD picture of the hadronic interactions. The identification of certain experimen- 
tally observed particles as glueballs will be a strong confirmation of QCD. Radiative 
J/\l/ decays have been regarded as the best hunting ground for glueballs IQ. Several 
candidates from experiments have existed for some time0. In order to compare the 
experimental findings with the theory, we need to find at least some characteristics of 
the glueball spectrum, or better yet to calculate some properties of the glueballs from 
QCD. Previous theoretical studies have concentrated on the estimation of the masses 
of the glueballs. The prediction of these glueball masses has long been attempted in 
lattice gauge Monte Carlo calculations . Recent progress using fuzzy loops 0, |^, H , 
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smeared loops [0, and inverse Dirac operators ^ to improve the signal to noise 
ratio have yielded quite consistent results on large volumes. These calculations show 
that the lowest scalar, tensor, and pseudoscalar glueballs lie in the mass range of 1 
- 2.5 GeV. In this letter, we shall report a lattice gauge Monte Carlo calculation of 
vacuum to glueball transition matrix elements which are needed for estimating the 
partial widths of J/\l/ radiative decays into glueballs. 

The glueball matrix elements we compute are of the form < 0|5'^G^^Gp^(x)|G' >, 
where G can be the 2"*"+, O"*"^ or ^ glueball and G^^ is the gluon field strength 
operator. On a lattice, the local operator g^G'^^,Gp^{x) can be constructed by small 
loops consisting of lattice links. One way to do the calculation might be, say for 
the 2++ glueball, to calculate the correlation function < O{t)O{0) >, where O = 
20i2 — Dsi — □23- Here, is a plaquette with edges in the and z/ directions. 
Then, for large Euclidean time, 

< O(t)O(0) >~ (< O|0|l >) V™^*, (1) 

where |1 > is the lowest glueball state with the quantum numbers associated with O. 
The matrix element can thus be extracted. Technically, however, we know from glue- 
ball mass calculations @]~[|lOl that small loops usually are not very good for creating 



glueball states. The states they create have small projections to the glueball states 
resulting in large noises in Monte Carlo calculations. In fact, it has become increas- 
ingly standard |^, |, |10[ in mass calculations to use the so-called "fuzzy operators" , 



which have extended structures, presumably reflecting the size of the glueball. On 
the other hand, it is difficult, if at all possible, to relate the fuzzy operators to the 
field strength operators numerically. Therefore, we adopt a method which is a hybrid 
of the two approaches. We use fuzzy operators to do a variational calculation first, 
finding the linear combination of different fuzzy operators that maximizes the over- 
lap with a glueball state in the calculation of the correlation function of the fuzzy 
operators < O{t)O{0) >. Then, we calculate the correlation functions of small loops 
with the fuzzy operator O 

< O{t)FF{0) >=J2<0\O\n><n\FF\0> e""""*, (2) 

n 

Here, FF is the lattice approximation of the G"'G"' operator. If we assume that the 
variational calculation gives a good overlap with the lowest glueball state as we let 
t » a, the lattice spacing, then we could fit the two correlations to the following 
expressions 

< O{t)O{0) >~< G|O|0 >2 e"™*, (3) 
< O{t)FF{0) >~< G|O|0 >< 0\FF\G > e"™*, (4) 

for different time t with three parameters, i.e. < G|O|0 >, < 0\FF\G >, and m. 
This way we could obtain the desired matrix element < 0\FF\G > and the glueball 
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mass m. Since these two correlation functions are measured on the same set of gauge 
configuraitons, we shall use the data-covariance matrix to fit them simulataneously 



11, 12 . More details to follow. 



Gluon configurations are generated by a pseudo-heat-bath algorithm []T3| . A mea- 



surement is made for every ten configurations generated. For each configuration to 
be measured, fuzzy links are constructed a la Michael & Teper0, as in Fig.l. The 
procedure produces links with length twice of the original links. It can be iterated 
several times to yield links of length 2^. In our calculations, 1 = 2. After this pro- 
cedure is done, the resultant links are no longer SU{?)) matrices. They do have the 
correct gauge transformation properties, however. In principle, these links can be 
used to form loop operators to create glueball states. In practice, however, we have 
found that they lead to huge statistical fluctuations. We use the following steps to 
convert the links back to SU{2)) matrices. Two conditions should be satisfied. First, if 
the fuzzy link happens to be an SU{?)) matrix, the procedure should keep the matrix 
un-altered. Second, the final matrix should have the correct gauge transformation 
properties. Let us consider a matrix M' and det M' ^ 0. We first divide it by its 
determinant M = M' / det M' . This obviously satisfy the two conditions. Now, any 
complex non-degenerate matrix can be factorized uniquely into a positive definite 
Hermitian part on the right and a unitary matrix on the left, M = UH [one can 
factorize the matrix as M = H'U' too resulting in a different set operators to create 
the glueball states). If under a gauge transformation M giMgl, by the unique- 
ness for the factorization, U ^ giUgl and H g2Hg2. Uniqueness also means that 
H = 1 and M = U if M is unitary. Therefore, when we convert M' to U, both 
conditions are satisfied. In actual computation, we first find the positive definite 
Hermitian matrix = M^M. We use the Jacobi method to diagonalize leading 
to H'^T = TA, where T is the orthonormalized unitary matrix of eigenvectors of 
and A is the diagonal matrix of eigenvalues. The inverse of H can be obtained from 
H^^ = TAr^/'^T^ . Finally, U = MH~^. This procedure is somewhat different from 
that used in gl]. 

To form the O operators for various spins, we use the U fuzzy links to form 
loops with shapes depicted in Fig.2. The same shapes were used in 0. Using the 
method outlined in |T5|, these loops and their rotations about different axes thereof are 
combined into operators belonging to the v4]^+(0++), E++(2++), T2++(2++), v4j"+(0-+) 
representations of the cubic group. The notations in the parentheses are the presumed 
quantum numbers for J^'-^ when the rotation invariance is restored. Only the real part 
of the trace is used because all the operators we consider here have positive charge 
conjugation, (a) and (b) in Fig. 2 contribute to all four kinds of representations, (c) 
does not contribute to the pseudoscalar. (d) and (e) only contribute to Af^{0^^) 
and E^~^{2~^^). These operators are averaged over spatial coordinates for each time 
slice to obtain the zero-momentum projections. 

The G^G" operators are constructed from the elementary(not fuzzy) links in Fig.3. 
Because of all the different shapes of the fuzzy operators and the yet-unknown varia- 
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tional coefficients, we do not know tfie polarization of tfie tensor states created by tfie 
fuzzy operators. We tlierefore need to calculate all the components of the FF tensor. 
For the pseudoscalar, we also need different components of the tensor in order for the 
FF pseudoscalar to have all the required symmetry properties. To get these, we first 
construct a lattice approximation of the field strength tensor. 

V = ^1 + ^2 + + U. (5) 

where Lj(i = 1,4) are products of link variables UUU^U'^ around the plaquettes in 
the plane as indicated in Fig. 3. 

= a'A,, + ^a\Dl + Dl)A,, + 0{a'). (6) 

Here, repeated indices are not summed over and A^^^, — igG^i,. a is the lattice 
spacing. Therefore 

Tr{F,^F,,) oc aYGl^G;, + 0{a% (7) 

This way, we can find all the components of the field strength tensor and all the 
Lorentz components of (gauge invariant) trace of products of the field strength tensor. 
The fuzzy operators are chosen according to the size and scale of the glueball. The 
operators for the field strength tensor, on the other hand, have sizes on the order of 
one lattice spacing. Once we choose the lattice scale such that the glueball sizes are 
much larger than the lattice spacing, then we expect that terms in higher powers of 
the lattice spacing a are relatively unimportant. This is the basis of our strategy to 
use large, fuzzy links to create glueballs and measure the field strength using small, 
elementary links. We shall use two {3 values to check the scaling of the results. 
A variational calculation which maximizes 

CjC, < O.(a)O,(0) > ^ < $|e-"^|$ > 

dCj < Oi{0)Oj{0) ><$!$> ^ ' 

is carried out to find the coefficients {q}. Here, |$ >= J2iCiOi\0 > and Oi are op- 
erators corresponding to different loops in Fig. 2. From the right hand side of the 
equation, we see that this is really not different from the usual variational calculation 
in quantum mechanics. Instead of taking the expectation value of the Hamiltonian 
itself and minimizing it, we take the expectation value of e~"^ and maximize it, which 
corresponds to finding the lowest eigenvalue of H. For the scalar case, vacuum expec- 
tation vahics arc first subtracted from the correlations, which corresponds to making 
the variational state |$ > orthogonal to the vacuum state. In actual calculation, this 
is equivalent to finding the eigenvector for the largest eigenvalue for the generalized 
eigenvalue problem, 

< Oi{a)Oj{0) >Cj = X< Oi{0)Oj{0) > cj. (9) 
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Once this is done, we proceed to evaluate the correlations < O{t)FF{0) > and 
< O{t)O{0) > using O = CiOi. Both the variational calculation and the fitting of 
Eqs. (0,|D are used to extract the desired matrix elements and glueball masses. 

For different spins, only the operators with same quantum numbers have non-zero 
matrix elements between the vacuum and the glueball states. Specifically, we would 
like to extract the following matrix elements for different quantum numbers 



Scalar 



Pseudoscalar 



s/a^ = <0\Tr{g'G,^G,^)\G> 

= 2 <0\Trg\E ■ E + B ■ B)\G > . (10) 



p/a' = <0\e^'^P'^Tr{g^G,,G,,)\G> 

= 4 < 0|Tr(/E ■ B)\G > . (11) 



Tensor 



te^u/a^ =< 0\e^,\G > (12) 

where e^^ = e^^, e^^ = 0, ef^^e^u = 1, and Q^^ = 2Trg^{G ^^Gua-l/ ^^^uG p^G p^) 
g'^T^jj where T^^, is the energy-momentum tensor. Here, 

\2 



t'/A g^a"^ = (< E,Ej > - < BiB^ >] 

+2< EiBj > {< EiBj > - < EjBi>). (13) 

and < * >=< 0| * \G >. 

The s, p, and t in the above equations are the lattice matrix elements obtained from 
eqs. (D and {§. 

Up to now, we have assumed the normalization < GIG >= 1. To go to the con- 
tinuum normalization < G, p|G, p' >= 2£'p(27r)^5^(p — p'), 

or < G,0|G , >= 2M{27rf6^{0) = 2MV, we need to multiply the above matrix ele- 
ments by V2MV. 

Calculations are performed for a 12^ and a 14^ lattices with (3 value equals to 5.9 
and 5.96 respectively. The lattice sizes are chosen so that the Liischer scale parameter 
z = m^++L |]T^ is about 10 where it is found [|1^] that E^^ and T2~^ start to become 



degenerate indicating the restoration of the rotational invariance. Furthermore, their 
physical sizes are about the same. From a cold start, we begin to make measurements 
after 4000 thermalization sweeps. The measurements are made for every 10 updates. 
For the 12^ lattice at /3 = 5.9, 5000 measurements have been performed. For the 14*^ 
lattice at /? = 5.96, 5800 measurements have been performed. Errors that enter the 
data-covariance matrix for the fit of eqs. (§) and (§) are estimated by binning 
the data into 20 sets and doing jackknife on them. The parameters chosen are in the 
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Table 1: Glueball Masses in unit of V^. ti/t2 denotes the range of time slices of eq. 
(3) / 6q. @ that are fitted. Also included are the per degree of freedom 





ti 


t2 


/3 = 5.9 


xVNdf 


ti 


t2 


/3 = 5.96 


xVNdf 




4,5 


1,2 


3.08(1) 


0.1 


1,2,3 


2 


3.42(1) 


7.9 


(0++) 


1,2,3,4,5 




3.29(9) 


0.14 


1,2,3 




3.40(3) 


0.04 


E++ 


2,3 


0,1 


4.44(1) 


19.2 


3.4 


0,1 


4.14(1) 


0.3 


(2++) 


1,2,3 




4.68(3) 


52.0 


2,3,4 




4.33(4) 


2.8 




1,2 


1,2 


4.83(1) 


190 


1,2,3 


1 


5.63(1) 


138 


(2++) 


1,2,3 




5.12(2) 


242 


1,2,3 




5.72(18) 


72 




1,2,3 


1 


6.10(3) 


49.7 


1,2,3 


1 


6.18(2) 


1.0 


(0-+) 


1,2,3 




6.05(3) 


40 


1,2,3 




6.17(2) 


0.8 



Table 2: Scalar matrix elements obtained from Eqs. ^ and ^. They are in the unit of 



(3 = 5.9 


P = 5.96 


ti 


t2 


s/a^ 




h 


t2 




xVNnF 


4,5 


2,3 


45.2(3) 


4.7 


1,2,3 


3 


51.1(4) 


7.94 


4,5 


1,2 


47.2(3) 


0.1 


1,2,3 


2 


53.0(1) 


7.9 



scaling(not necessarily asymptotic scaling) regime^]. We present the results in units 
of the string tension K, which is determined from fitting the correlations of Polyakov 
loops to the exponential Ae~^^* for the time interval t = 2 to 4 for /5 = 5.9 and t = 
2 to 5 for /5 = 5.96. The scale determined from the string tension indicates that the 
lattice spacing decreases by ~ 19% from /? = 5.9 to 5.96. This makes the physical 
sizes of the 12^ and 14^ lattices about the same. The glueball masses determined from 
both the joint fit of the correlations of eqs. (^ and @ and the fit of the correlation 
of fuzzy loops alone in eq. (^ are listed in Table 1. They are given in units of 
\/K. They agree with each other and with the higher statistics calculations in this 
f3 range ||^ (i.e. /? = 5.9, 6.0) for each mass determined from the corresponding time 
slices. However, for the T2"^{2~^^) case, we could not find a fit with small x^- This 
presumably reflects the fact that our lattice is still not large enough for the tensor 
glueball, since it is known that the tensor glueball is about 4 times larger than the 



scalar glueball |T6|, |17 



The glueball matrix elements with several fits to eqs. (^ and (§) are list in Tables 
2 through 5. They are in the unit of K^^^. 

These numbers can be converted into physical units using K ^ {0A2GeV)'^{i.e., 
multiplying the numbers by K^^"^). We see from these results that the scalar glueball 
matrix elements are better calculated, much like the calculation of the glueball masses. 
The variation from f3 = 5.9 to 5.96 is about 10%. This is quite reasonable as the 
matrix element scales like K^^'^. The rest of the matrix elements differ more. Again, 
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Table 3: Tensor matrix elements for the E~^~^ representation in units of K^^'^. 



(3 = 5.9 


P = 5.96 


h 




t/a* 


xVNdf 


h 




t/a^ 


xVNdf 


2,3 


0,1 


7.42(6) 


19.2 


3,4 


0,1 


11.0(1) 


0.3 










2,3,4 


2 


8.64(14) 


4.5 



Table 4: Tensor matrix elements for the T^~^ representation in units of K^^'^. 



(3 = 5.9 


p = 5.96 


h 




t/a^ 


xVNnF 


h 




t/a^ 




1,2 


1,2 


6.38(8) 


190 


1,2,3 


1 


10.1(1) 


138 



the case for Tj*""^ is the worst with no fit of small x^, enough though the values of 
the matrix element are not that dissimilar to that of Both E^^ and Ai^{0^^) 

have fits with reasonably small x^- 

To come to grips with the continuum matrix elements, we need to consider the 
finite lattice renormalization of the lattice operators. The finite renormalization of the 



lattice operators considered here have been calculated perturbatively |]T8|, |1^, ^ 
We can define the renormalization of the operators as follows, 

O = ZOia-'^ (14) 

where the 0^ is the lattice-regularized version of the continuum operator O and d is 
the mass dimension of O. For the case of the scalar operator g^G^yG^^ which coincides 
with the renormalization group invariant trace anomaly ^^^^G^uGfj,^ to lowest order 



in g, the renormalization constant Z is found [0 to be 1.06 at one- loop order. For 
the pseudoscalar lattice operator which is the topological charge in the limit of zero 
lattice spacing, the finite renormalization Z is found to be 2.5 with the two- loop 
corrections estimated by dominant tadpole graphs [0. For the energy- momentum 
tensor operator, the lattice perturbation to one loop gives the Z to be 1.84 for j3 = 5.9 



21 



When these renormalization constants are taken into account, the scalar matrix 
element in eq.(0) is about 50^3/2 or about 3.7GeV^ for K = {QA2GeVf. Based 
on the scaling properties of QCD and trace anomaly, both the QCD sum rule |^ 
and the soft meson theorem lead to the estimate that relates the scalar glueball 



Table 5: Pseudoscalar matrix elements in units of K^^'^. 



(3 = 5.9 


(3 = 5.96 


ti 


t2 




xVNdf 


ti 


t2 




xVNdf 


1,2,3 


1 


19.2(2) 


50 


2,3 


1,2 


21.9(2) 


210 


2,3 


1,2 


17.4(2) 


0.05 


1,2,3 


1 


23.2(2) 


1.0 
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matrix element to the gluon condensate, 




< Q\Tr{g'G^,G^,)\G >= lQn\ -^mc, (15) 



where Go =< Ol^G^G'^.^lO > is the gluon condensate, b = 11/3N^ - 2/3Nf and 
the scalar glueball mass. Taking Nf = for the pure gauge case, Gq = 0.012GeV^ 
from the ITEP value [^, and mc in the range of 1 to 1.7 GeV, this matrix element 



is then estimated to be 3.7 ~ 6.3 GeV. This agrees with our lattice calculation to 
within a factor of two for the extreme cases. For the pseudoscalar case, the lattice 
renormalization Zp = 2.5 makes the glueball matrix element in eq. (^) to be about 
or 3.3 ~ A.SGeV^. It has been proposed that there is an approximate 
chiral symmetry between the scalar and pseudoscalr glueballs A sum rule is 

derived from an effective action which relates the topological susceptibility x the 
pure gauge case to the gluon condensate Go 

(16) 



where ^ denotes the degree of chiral symmetry. From the Witten-Veneziano mass 
formula for the r]' mass which is related to x the large Nc analysis, this sum rule 
predicts ^ to be 0.7. From our lattice results, ^ can be obtained from the ratio of the 
pseudoscalar to scalar matrix elements 



- 2z:ij^^- ^^^^ 

When we take the ratio between the pseudoscalar and scalar matrix elements from 
Tables (2) and (5), C,l truns out to range from 0.5 to 0.6. This is quite close to the 
prediction from the sum rule in eq. (0). In the Witten-Veneziano mass formula for 
the rj' mass [^], the glueball contribution to the topological susceptibility is neglected. 



The glueball contribution to the topological susceptibility is 

_ <0|rr(^^GG)|G> , 

Xglueball — [ -Tjr^ J l-LoJ 

Using our lattice result and the mass of the pseudoscalar glueball at 1.4 GeV, the 
pseudoscalar glueball turns out to give an appreciable 9 % contribution to the topo- 
logical susceptibility. With the perturbative lattice renormalization Zt-, our lattice 
result for the tensor glueball matrix element is ~ lAK^^"^ or ~ 1.0 GeV^^. On the 



other hand, the prediction from the tensor dominance model |2^ and QCD sum rule 
is about 0.16 ~ 0.35 Gey^ for the tensor glueball mass ranging from 1.7 to 2.2 
GeV. Hence, the lattice result is about 3 to 6 times larger. Since the tensor dom- 
inance model gives an reasonable prediction of 6'(/2(1720)) production rate in J/\& 



radiative decay we suspect that the lattice result we obtain for the tensor 
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case is overestimated. Recent glueball wavefunction study |T6|, |T^] reveals that the 
size of the tensor glueball is about 4 times larger than those of the scalar glueball 



and the pion Hence, we expect that the tensor glueball is more susceptible to 
the finite size effect. Considering that our choice of the lattice size should be large 
enough for the scalar and pseudoscalar glueballs but may not be enough for the ten- 
sor, the squeezed tensor glueball will lead to an overestimated tensor matrix element. 
In the constituent glue picture, the matrix element we calculate can be viewed as the 
wavefunction at the origin which scales as the inverse 3/2 power of the glueball size. 
We think this is the likely source for the large tensor matrix element. 

Finally, we would like to discuss the prediction of the glueball production rates in 
J/\E' radiative decays. While the details of the calculation will be presented elsewhere 
pl| , we would like to mention some of the main features of the findings. It is shown 
that the partial widths of the glueball production in J/^> radiative decays can be 
related to the vacuum to glueball transition matrix elements calculated in this paper 



27| , |31[1 . This is based on the approximation that the charm quark is heavy so 
that c and c annihilate at a point. The recent analysis of the vrvr and KK decays 
shows that a large component of spin zero is observed in the ^(1720) region [p2| . 



If this is confirmed, it would be a very good candidate for the scalar glueball [pO 
The calculated scalar glueball matrix element in eq. (p!0| ) predicts a branching ratio 
of 5 X 10~^ in J/^> radiative decay which is a only a few times larger than the 
the experimental value from the observed decays in vrvr and KK channels |3^. The 
unobserved decay modes will make the agreement better. The calculated pseudoscalar 
matrix element gives a prediction of 4 x 10~^ for the branching ratio of the 0~ glueball 
candidate ?7(1440) in J /'^ radiative decay. This agrees quite well with the observed 
B.R. of ~ 3.5 X 10-3 form the decay modes of KKTi,KK,aQ{dmy, and pV° [§3- 
On the other hand, the predicted branching ratio for the tensor case, assuming that 
the 2"*" component in ^(1720) is the tensor glueball, is one to two orders of magnitude 
larger than the experimental one. This could be due to the fact that the tensor 
glueball matrix element is overestimated as alluded to earlier. Or perhaps the tensor 
glueball lies higher in mass (e.g. > 2 GeV) leading to a smaller predicted B.R. 

To conclude, we have calculated the glueball to vacuum transition matrix elements 
in quenched lattice QCD. Notwithstanding the fact that our present study does not 
have high statistics and the Euclidean time separation is limited, the variational 
method for the hybrid correlation functions of the fuzzy and plaquette operators does 
seem to yield consistent results for the glueball masses and matrix elements. The 
matrix elements obtained for the scalar and the pseudoscalar glueballs seem to agree 
reasonably with those obtained phenomenologically and the prediction for glueball 
production rates in J/'^ radiative decays are in line with the experimental results 
based on the observed decay modes. The tensor case could be an exception. It yields 
too large a branching ratio for the tensor glueball in the 9 region in J/\E' radiative 
decay and is considerably larger than those from the tensor dominance model and 
the QCD sum rule. We speculate that it is due to the large finite size effect which is 
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particularly acute for the tensor due its large size |[T^, |T^. To verify this, one needs 
to conduct a high statistics calculation on a larger lattice which should restore the 
rotational invariance for the E^~^ and T2"^ and allow a larger separation in t in order 
to have finite size effect under control. 
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Figure Captions 

Fig. 1 The making of fuzzy links a la Michael and Teper. 

Fig. 2 Shapes of different loops used in the variaitonal calculation. 

Fig. 3 Operators for the local field strength. 
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